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Abstract 

We discuss a model based on a field-induced mixture of two odd-parity irre- 
ducible representations to explain the unusual features of H C 2(T) in the heavy 
fermion compound UBei3. We compare its predictions with recent pressure 
measurements as well as with the most prominent theoretical models which 
have been proposed up to now. 
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I. INTRODUCTION 



Heavy-fermion compounds are well-known candidates for a search and investigation of un- 
conventional superconductivity. In these three dimensional systems, the symmetry analysis 
of the possible superconducting state, depending on the crystalline lattice of the compound 
ffl, sets the frame for the identification of the unconventional superconducting phase. This 
analysis gives a complete list of possible superconducting phases together with their proper- 
ties determined by symmetry, including the type and order of nodes in the superconducting 
gap. In spite of a favorable theoretical situation and numerous experimental work, there 
doesn't exist yet a firm and widely accepted identification of these superconducting phases 
(comparable to the results on superfluid 3 He) in any heavy fermion superconductor. 

In this paper, we examine the upper critical field H c2 of the heavy fermion superconductor 
UBei3 as a probe of the symmetry of its superconducting phase. The upper critical field is 
usually not very sensitive to gap nodes (see for example the quantitative studies for various 
scenarios in UPd2Al 3 0). But it is sensitive to the spin state of the Cooper pairs or more 
generally to the parity of the order parameter, because in heavy fermion systems, the orbital 
limitation is so large that may be governed by the paramagnetic limitation. 

UBei3 is a cubic compound with a T c of order IK, and it is a non-magnetic superconduc- 
tor. Since the first measurements of its upper critical field [Q, it is known that it presents 
two very unusual and intriguing features. First, H C 2(T) has a strong negative curvature 
close to T c which changes sign at intermediate fields. Then, taking account of a realistic 
value of the conduction electrons gyromagnetic g — factor, the paramagnetic limit at T=0 
is exceeded several times, while the strong negative curvature close to T c shows that there 
exist a pronounced effect of the paramagnetic limitation. 

Numerous explanations have been proposed since the first precise measurements of 
H C 2(T) in this compound ||. Some have relied on additional hypothetic magnetic phase 
transitions (4|, or on the field dependence of the normal state properties But none 

of these phenomenological interpretations have found a firm basis in other measurements 
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or theoretical developments. Another hypothesis, much closer to the point of view adopted 
here, relies on two different superconducting order parameters with a weak Pauli limitation 
||. It has not been carried out quantitatively, but experimental support has been seeked 
through the detection of a possible second phase transition already in zero field below the 
main superconducting transition To our opinion, such a second phase transition is 

not supported by the data, which only show a weak and smeared maximum in the specific 
heat or minimum in the thermal expansion. The model that we propose here involves a 
mixture of two different irreducible representations, but does not rely on nor predict such a 
second phase transition. 

At present, the only competing quantitative explanation relies on a simple strong- 
coupling model |PH| . It has been also successful in describing the evolution of both features 



up to pressures of 20 kbars [|1^]: the complete temperature and pressure dependence of H c2 
comes out from a straightforward strong-coupling calculation with a single even parity state. 
The conflict with the paramagnetic limit at T=0 is resolved by its enhancement due both 
to direct strong-coupling effects (increase of the ratio A/T c ) and to the (parameter free) 
inclusion of the formation of a spatially modulated superconducting state (FFLO), induced 
by the dominance of the paramagnetic limitation. Very good agreement with the data of 
Ref. [1^] is provided in the whole pressure range, for g being close to its free electron value, 



by fitting the strong coupling constant A. The pressure dependence of A agrees with that 
of the effective masses (as indicated by the Sommerfeld coefficient or the slope of H c2 at 
T c ), but it also turns out to be exceptionaly large: A ~ 15 at p=0. For all other known 
superconductors where the strong coupling regime due to electron-phonon interaction is well 
characterized, A does not exceed 5. 

The model in this paper is based on a field-induced mixture of two odd-parity irreducible 
representations of the symmetry group {Oh) of the normal phase of UBei3. This possibility 
was already evoked in the literature [p!3| , p^| , but it has not been worked out to the extent 
that it could justify the particular choice of the order parameters and could be compared 
directly with the experimental data. Here, the choice of the two representations is unique, 



thus providing a firm identification scheme for the superconducting state of UBei3. Within 
this scheme, the dominant component of the order parameter is analogous to the B-phase of 
the superfluid 3 He. We make a quantitative comparison of the theoretical predictions with 



the data of Ref. JT2| and propose an experimental check of the suggested scheme. 

As we shall see, our choice of order parameters does not allow for a consistent interpre- 
tation of all data in UBei3: in particular, it does not predict any nodes of the gap, whereas 
both old [16| and new fl7| thermodynamic experiments hint at the presence of point nodes. 



We consider this approach as a first step which may help to clarify at least one aspect of the 
properties of UBei3, bearing in mind that the complexity of the physics of heavy fermion 
systems rarely allows for fully satisfying explanations. . . 

II. CHOICE OF THE REPRESENTATIONS. 

Rotational symmetry of the crystalline lattice of UBei 3 is described by the Oh group. 



Following the list of Ref. JT8[ for strong spin orbit coupling, we have to consider 10 irreducible 
representations: Ai 9iU ] A 2 g jU ] E gu ; Fi g>u ] F 2g)U . The subscripts g and u denote correspond- 
ingly even or odd symmetry of the particular representation with respect to the inversion. 
The even representations are the analogs of the spin singlet pairing and the odd - of the spin 
triplet. The order parameter for an even representation is a scalar and for the odd - a vector 
function of the direction in the momentum space d(k). The capital letters denote A - one- 
dimensional, E - two-dimensional and F - three-dimensional representations. To determine 
which of the representations is realized as the order parameter in a particular compound, 
one has to know the projections of the scattering amplitude V a p\^{k^ k') of quasiparticles in 
the normal phase on the basis functions of these representations. The scattering amplitude 
is then represented in the form: 

z r 9 i=i 

+~ E Vr, Ei^ikUM^fik).^) (1) 
z r u j=i 



The spin matrices are defined as g a p = (iaa y ) a/ 3, g a p = io y afi where CTaB^aB^aBi are Pauli 
matrices. The sums are taken over all even and odd representations respectively. The 
greatest of the positive amplitudes Vr (the dominant one) determines the temperature of 
the real superconducting transition T c0 as prescribed by the BCS-theory 

where In 7 is Euler's constant, er a cut-off parameter and N(0) the density of states at the 
Fermi level. All other amplitudes are subdominant. The transformation properties of the 
order parameter not too far from T c0 are described by the representation T . A magnetic 
field (H) changes the symmetry of the system and the classification of the representations. 
If only one, dominating term Vr is kept in the sum (|1|) and if T is a multidimensional 
representation, the magnetic field splits it. Several branches H C 2(T) start in that case in the 



plane (T,H) from the point (T c o,0). According to Ref. |Tj|, if the direction of the field H 
coincides with the symmetry axis of the lattice, the different branches are classified by two 
quantum numbers. The first is N=m+n, where m is the projection of angular momentum 
on the direction of field and n is the number of the Landau level, which describes the spatial 
dependence of the order parameter for the nucleating superconducting region. The second 
quantum number is the parity a = ±1 with respect to reflection in the plane perpendicular 
to the direction of H. The branch with the highest H for a given T determines the true 

The dependence H c2 (T) for p-wave pairing has been analyzed in detail by Scharnberg 



and Klemm |19| . No change of sign of the curvature was found, which is natural in a problem 
with only one parameter Vp . But if all terms are retained in the sum ([]]), the magnetic 
field has an additional effect: it can mix basis functions having the same quantum number 



1 Strictly speaking even for conventional s-wave pairing, there are many branches H C 2(T) starting 
from T c and corresponding to different n, but one knows that the highest branch corresponds to 
n=0. 



N, which originally belonged to different representations. Substantial change of the original 
basis functions can take place for fields of order ^bH ~ T c . A change of the basis functions 
will change the projections V a pxn(k, k') on these functions, including the projection which 
determines the superconducting transition. A substantial admixture of the subdominant 
coefficients Vp may therefore also take place on energy scales ~ T c , imitating an extreme 
strong coupling effect. As regard the field induced change of the interaction potential itself, 
it is of the order of [XbH/cf an d can be neglected. 

The minimum model which can include these effects must contain two representations 
well mixed by the magnetic field. An advantage of the minimum model is that it re- 
mains tractable and contains only one additional parameter. Whether or not the two- 
representations model is sufficiently accurate for a description of the magnetic properties of 
a particular superconductor depends on the values of coefficients in the sum ([I]) for this ma- 
terial, and on the temperature interval in which one expects to reproduce these properties. 
In what follows, we apply the two-representation model to reproduce the unusual temper- 
ature dependence of H c2 in UBei3. We use the above mentioned characteristic features of 
this dependence as a guidance for the choice of two dominating representations. 

The critical magnetic field can exceed the paramagnetic limit if the order parameter 
in the high field region is predominantly of the odd-parity type (i.e. a vector-function 
d(k)), and if the projection of d(k) on the magnetic field is small or absent. To provide a 
good mixing of the two participating representations by the magnetic field we assume that 
both representations are of the odd-parity type. The pronounced paramagnetic effect in 
the low field region indicates that in that region the order parameter is dominated by a 
one-dimensional representation. In that case it is not possible to eliminate completly one 
projection of d(k). Such an elimination becomes possible in strong fields when the admixture 
of the second representation is appreciable. An inspection of the representation table for 
the Oh group with the corresponding basis functions |2D| suggests the following scheme: the 
superconducting phase which appears at H=0 and T=T c0 belongs to the A\ u representation 
with the basis function 
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$o(k) = xk x + yk y + zk z 



(3) 



The unit vectors x, y, z are directed along three mutually perpendicular four- fold axes of 
the cube. The order parameter of the form (|3|) corresponds to a state of spin S=l, orbital 
momentum L=l and total angular momentum J=0. This is the most symmetric odd-parity 
state, which is analogous to the B-phase of superffuid 3 He. The superconducting gap for that 
state has the full symmetry of the Oh group and it does not have nodes required by symmetry. 
We assume in what follows that the magnetic field is oriented along the z-direction. Only 
one branch of //^(T) starts from T c o (ignoring higher Landau levels) and this branch is 
characterized by the quantum numbers N = and a = +1. The second representation is 
E u . For the present discussion it is convenient to choose the basis functions of E u in the 
following form: 



is zero and the paramagnetic limitation is absent for such an order parameter. The elimina- 
tion of the z-component would not be possible for the Am representation, since as a function 
of k it is orthogonal to z-components of the basis functions of all other odd representations. 
So, we assume that the pairing potential contains contributions of only two representations 
A lu and E u . 



where e 




with m = ±2. One can see immediately that the z-component of the combination 



= v / 2(l + e)* + ^ 



(6) 



VWMO = -v (Mk).g^)(%(k).^) 



(7) 



8=1,2 
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with two independent coupling constants Vq and V\ or, formally, two BCS transition tem- 
peratures T c0 and T cl . Then the order parameter in the vicinity of the transition line H c2 (T) 
is a linear combination of the basis functions ^o, ^i, ^2- 

d(k, R) = A (R)$ (k) + ]T &s(R)$s(k) (8) 

8=1,2 

where Aj(R) j=0,l,2 are functions of the coordinates. 

III. CRITICAL FIELD H C2 (T) 

The dependence of if c2 (T) for our model can be found as in Ref. |19j. As a starting point 
we use Eq.(3) of Ref. JHJ. With the interaction potential given by Eq. (^) and the order 
parameter of Eq. (Q), the equations for the functions Aj(R) are: 

1 2 r dVt' - r°° 

e~ sLop .{l - [1 - cos{2sgH)}nn tr }A j/ (R)^ f {k'), (9) 

where 

L op = 2\u n \ - isgn(u; n )v F (k F ).(iV + 2eA(R)) (10) 

j = 0,1,2; A(R) is the vector potential for the magnetic field H , g is an effective gyromagnetic 
ratio. It need not be equal to its free electron value g = 2: besides the usual renormalization 
due to spin orbit coupling, it includes here also possible Fermi-liquid corrections. In what 
follows, g will be used as a fitting parameter. The unit vector h is parallel to the H direction 
(thereafter also called z-direction). This is a linear system of differential equations of infinite 
order. Its solutions are given by the eigenfunctions f n (R) for the Landau levels of an electron 
under magnetic field: 

A 1 (R) = J2v 1 nUR) (11) 

n 

Now we use the fact that the combination (^) corresponds to the quantum number N=0. 
For the functions \l/o an d ^i m=0, while for \& 2 m = ±2. This selects n=0 for j=0,l, and 
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n=2 for j=2. Since the spatial dependence of f n (R) is known, we suppress the second index 



in the notation r\j n and after transformations following that of Ref. ||19|| , we arrive at a linear 
algebraic system for the amplitudes r\f 



(F 00 + P - In Vh)7] + V2(l + e)(F i - P)r]i 
1-e 2 



+ 



V2 



-F 02 T]2 = 0, 



V2(l + e 2 )(F 01 - P)Tjo + (F 00 - F 01 + 2P 
- In Vh - In q)rji + ^— - -F02772 = 0, 



1-e 



1-e 2 



Fo2?7i 



+ (F 00 + F i - In - In q)r) 2 = 0. 



(12) 



(13) 



(14) 



The equations are written in the dimensionless units of Ref. [|19j: h — H — 



2ttT c , 



T 

TcO '■ 



, $o — t = q = For the calculation of angular averages, we assume a 



Is 



spherically symmetric Fermi surface so that Vp is parallel to kp. Then the coefficients F Q0 , 
Fqi, Fq2 and P are given by the following expressions: 



oo 



'01 



02 



00 



.Vh, 



In 



Vh ( pt s 

tanh — = 

t \2y/h. 



x 2 )e~ 



'01 



.Vh, 



Vh 



pdp 



In 



o 



— tanh ( — = 
* \2Vh, 



o 2 , 



pdp 



t 



t 



P 2 dp 

iY2 \Vh) 4V2hJo sinh(^ 

(Jjl-x^e^-^dx 

Vh Vh Jo smh(^) 



l^e-Td- 2 )^) dp 



(15) 



(16) 



(17) 



(18) 



with A = f^rf"- H C 2(T) is found from the condition of compatibility of Eqns. (|T^|)-(|rj|). 
With the shorthand notations 
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G = F q — In yh — |Q, where Q = In q, 



this condition has the form: 



(G + ^)(G-|) 2 -F 2 1 (3G + 2F 01 ) + 

3P[(G + Foi) 2 - ^] + ^ 2 2 (2F 01 - 3P - G) = (19) 



which gives an implicit equation for h. 



IV. DISCUSSION 

Analytic solution of Eq. ( |19"D is possible only in some limiting cases (cf. Appendix). In 
order to compare the predictions of Eq. ( |T9"D to the data, we have performed straightforward 
numerical calculations. Fig. [I] shows H C 2(T) for three different irreducible representations 
together with the data on UBei3. A pure A\ u representation has a paramagnetic limitation 
3 times higher than a pure (even) Ai g , leading to a good fit of the low field part of H c2 (T) 
for a reasonnable value of the effective g — factor: g = 1.2. But the admixture of the E u 
representation is essential to reproduce the upturn of H c2 (T) below T c /2. 

The surprise has been that this large influence of the admixture has been found for a 
parameter T cl /T c0 of only 0.12. This comes from the fact that the admixture of the repre- 
sentations along the H c2 {T) curve is controled by the applied field, and it starts much above 
T cl . In any case, the corresponding admixture of E u components in the order parameter 
(defined in Eq. (|TTD) remains less than 10% down to T=0: see Fig. [|, where the coefficients 
7]i and r)2 of definition (|8|)-(|TT|) have been reported normalized to J2Vi = 1- So one is still 
far from the limit of complete suppression of the paramagnetic limitation (see Appendix). 
One can note on Fig. ^ that most of the admixture of the E u representation comes from the 

function (Eq. (|J)), which compensates the m z = component of the A\ u representation 
with a Landau level n = instead of n = 2 for the ^2 function (Eq. (^)). 

A complete comparison of the best fits obtained in our model with the results under 
pressure of Ref. JL2] are presented on Fig. |3|. Three parameters are used in the process of 
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fitting: vp, which is directly found from the slope (— ^ a J T T , the effective g — factor, 
which is determined by the curvature of the low- field part of the data, and the ratio |fj, 
which is controled by the high-field region. 

The values of the fitting parameters for different pressures are presented on Fig. 
They suggest several comments. First of all, we note that the two- representations (Ai u , E u ) 
model can reproduce the unusual features of H c2 (T) in UBei3. Good quantitative agreement 
is obtained everywhere except in high fields. At least two reasons can be invoked for the 
observed deviations: the effect of the omitted representations in equation (0) is expected to 
become stronger in the low temperature range, and the strong coupling effects, known to 
exist in this compound ( 0,^]), have been neglected in our calculations and would also 
reinforce H c2 (T) in this temperature range fl2f . So these deviations can be ascribed to some 
oversimplifications of our model, without fundamentally questionning its validity. 

The fitting parameters have realistic values. The pressure dependence of vf is in agree- 
ment with that of the Sommerfeld coefficient (deduced from Ref. |21]). This is not a surprise, 



as this feature is almost model independent and it has been already pointed out in Ref. [12 



The effective g — factor increases slowly from 1.2 to 1.5 and is not too far from the free 
electron value. In the absence of any other probe of g, no additionnal cross-checking seems 
possible. 

The more characteristic parameter of the model is the ratio T c \/T c q. Pronounced effect 
of the subdominant interaction via admixture of the second representation is observed even 
when the temperature of the subdominant transition T c \ is much smaller than the real 
transition temperature T c0 . At ambient pressure T c i/T c0 is of the order of 0.1, whereas 
it grows linearly with pressure in the interval < P < 20kbar, up to a value of 0.7. 
Extrapolation of this dependence in a region of higher pressures predicts a crossover of the 
two temperatures at P « 30kbar. For P > 30kbar, T cl would be larger than T c0 and the 
transition in zero magnetic field will take place in one of the three possible phases belonging 
to the E u representation |L8| . An experimental observation of this crossover would be a good 
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check of the proposed two-representation scheme. The details of the crossover (if observed) 
could give support to the proposed choice of the two representations. 

The observed qualitative changes (suppression of the upturn) of the curve H c2 (T) with 
increasing pressure are easily understood within the proposed scheme. When T c \ approaches 
T c0 , two representations can be considered as one three-dimensional representation, in which 
case the paramagnetic limitation can already be suppressed in low fields by a suitable com- 
bination of the basis functions, and no upturn will appear. A possibility to fit better the 
experimental curve when T& (and Vi) is growing is also natural, since in that case the role 
of the omitted representation is getting relatively smaller. 

Eventually, let us note the difference between our model and the two transition proposal 
of ref. ||. In the latter one, it is proposed that a second transition is present below m 0.6T C , 
and that the upturn comes from the apparition of a new phase with a weaker Pauli limi- 
tation. The recent experimental effort reported in [|1(J strives to detect such a transition 
already in zero field. In our model, the second phase (the E u representation) does not appear 
in zero field-zero pressure. It is only the applied field which introduces the mixture of the 
representations at finite temperatures (see Fig. ^). A new phase transition is avoided pre- 
cisely because the field allows such a mixture without additionnal symmetry breaking. The 
present scenario, in addition to providing quantitative predictions, also has the advantage 
of giving a realistic account of the smoothness of all features observed below T c0 in UBei3. 

V. CONCLUSION 

The proposed scheme satisfactory explains the temperature dependence of H c2 (T) in 
UBei3, but it is in conflict with the observed power-law temperature dependencies of the 
specific heat, London penetration depth and of the longitudinal relaxation time in NMR- 
experiments [[HJ which reveal the existence of gap nodes: with the A\ u order parameter, 
the nodes in the superconducting gap can only be accidental. The argument given above 
for the choice of the representations indeed leaves no other possibility than this A\ u — E u 
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mixture among odd representations. There remains a possibility of a mixture of even and 
odd representations [p~^JT5f| . The choice of possible pairs of competing representations in that 



case is much bigger. We have not analyzed that possibility systematically though one would 
expect that in that case, because of a weaker coupling between the representations, the 
change in the curvature would be sharper. A promising route could also be the inclusion of 
strong coupling effects, which have recently received microscopic experimental confirmation 
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APPENDIX: 

We consider now some limiting cases where the solution of Eq. ( p~9|) with respect to h 
can be found analytically. This analysis is useful both for a clarification of the underlying 
physical picture as well as for a check of the numerical calculations. Let us first consider 
the limit T — > 0. The limiting values of the functions of the variable t/y/h F 00 , Fox, F 02 for 
are: 

Foo( 0) = r M ^)^/ 1 (l-^)e-^) lfe 
Jo 2 2 Jo 

= ln(^=) (Al) 
P\Pdp 



*■„«» = jf ln(f, 4 

\l-x 2 )(l-?,x 2 )e-^ 1 ' x2) dx 
1 

~ ~6 

F 02 (0) = —7^1 Pdp J Q (1 " J)e-W-* ] dx 



(A2) 
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P(4=, \Vh) in a limit t/Vh -> remains a function of A\//^ 
P = P(0, Av 7 ^) 

= _ P ^£[1 _ cos(pA^)] ^V^ 1 "* 2 ^ (A4) 
jo p Jo 

For a further simplification we consider the limit \Vh — > oo as being relevant to the actual 
situation in UBei3. To evaluate the asymptotics of the integral in Eq. ( |A4j) in that limit, 
let us split the interval of integration in two: (0,po) an d (Po, oo) where po is chosen to meet 
the following condition: 

AV^ < P ° < 1 (A5) 

In the first interval, one can then assume e~~^~ x ) ~ 1, and the contribution Jx of this 
interval is evaluated straightforwardly: 

7 1 rfp cos(pAv ^ }] = 1 l n ( 7 Av^po), (A6) 

o JO p 6 

where In 7 = C is Euler's constant. When integrating over the second interval (po, 00), one 
can drop the oscillating term cos(p\y/h), yielding for the contribution J 2 of this interval: 

h = Jo ^ C f e~«, where u = (1 - x 2 )^ 

Integrating over u by parts and taking into account that 

/ In ue- u du = -C = - In 7 (A7) 

JO 

we arrive at the following contribution of 1%: 

J 2 = -iln(p 0v /7) + ^ (A8) 



The sum of ( |A6[ ) and ( |A5| ) gives the principal order terms in the asymptotics of P for 
\\/h — » 00: 
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(A9) 



Given the limiting values Eqs. ( |A1| ) - ( |A3| ) and the asymptotic (|X9|), one can find a limiting 
value of H C 2 at T=0 and Av^ — > oo. According to ( |A9| ) P — > oo when Av^ — > oo. Collecting 
in Eq. (^) the terms proportional to P and setting them to zero we arrive at the following 
equation: 



(G + F 01 ) 2 = ^ + ^ 2 2 



(A10) 



Using here the limiting values Eqs. (|A1|) - ( |A3| ) we solve this equation with respect to h. 
The largest of the two roots gives H c2 : 

1 



4 7 g4/3 



cxp 



5 + 2,^ + 1 
3 V 9 8 



(All) 



For this value of h one obtains from Eqns. ([12]) - flH ) the limiting values of the ratios of the 
coefficients in the definition (O) : 



771 1 + e 2 



(A12) 



Vo V2 

This ratio corresponds to a complete elimination of the z-projection of the order parameter 
(cf. Eq. ©). 

Another ratio is: 



1-e 



Q Q* 1 
3 + VT + 8 



V2 3 

The order parameter obtained with these coefficients has the form: 



(A13) 



d(k) = VSyky 



(A14) 



In the opposite limit when T — > T c and H C 2 -> we can expand the functions entering 
Eq. flT9"D in powers of h and r = 1 — t. These functions are defined as double integrals. For 
integration over p the convergence is provided by the exponential factor e ^ in a region 
Vh/t < 1. It means that all other functions can be expanded in powers of p 2 which, 
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after integration over p, gives an expansion over h/t 2 . We keep only terms which are neces- 
sary for finding r(h) with an accuracy up to h 2 : 



F 00 - InVh = -hit - a£ + &£, F 01 = -f^, 

with a = ^C(3), & = §C(5), fi = ^C(3), h = ji^Ci 3 ) where C(*) is Riemann's zeta 
function. Substitution of this expansion in Eq. ( |19|) gives 

r(/i) = a/i + (L 2 -b + a\ 2 - i[2/ 2 + ^/l])/* 2 (A15) 

The term linear on h comes from Fqq and is determined entirely by the A\ u representation, 
but the h 2 term is influenced by the E u representation via the functions Fqi and Fq2 (or the 
coefficients fx and Z^)- This influence is getting stronger when Q decreases, i.e. when the 
two transition temperatures are coming closer. The term proportional to 1/Q is definitely 
negative: it means that for sufficiently small Q, the dependence of H c2 (T) will have positive 
curvature starting from T c . 
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FIGURES 

FIG. 1. Upper critical field of UBei3 at zero pressure (from [0]), together with three calcula- 
tions of H C 2{T) for three different irreducible representations: A\ g , A\ u , and a mixture Ai u — E u . 
The same values of the effective g — factor (adjusted for the odd parity representations) and the 
Fermi velocity have been used for the three fits. 

FIG. 2. Temperature dependence of the normalized components r/i and 772 of the respective 
basis functions ^1 and ^2 of the E u representation, along the H C 2(T) fit of Fig [|. Note that they 
both start to grow much above T c % = 0.12T c q. 



FIG. 3. Upper critical field of UBei3 under pressure [11], and best fit of our A\ u — E u model 
(full lines) for each pressure: the main features are well reproduced, and deviations appear only at 
low temperatures (see discussion in the main text). 

FIG. 4. Pressure (p) dependence of the three parameters of the fits of Fig. ||. The comparison 
of the p-dependence of the Fermi velocity vf with the Sommerfeld coefficient j(p) of ref. |^] is 
gratifying. Extrapolation of the linear increase of If 1 predicts that the E u representation will 

J- cO 

appear first in zero field above p « 30kbar. 
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